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Abstract 

We present the results of the modeUng of proton translocation in finite H-bonded chains in the 
framework of two-stage proton transport model. We explore the influence of reorientation motion 
of protons, as well as the effect of electric field and proton correlations on system dynamics. An 
increase of the reorientation energy results in the transition of proton charge from the surrounding 
to the inner water molecules in the chain. Proton migration along the chain in an external electric 
field has a step-like character, proceeding with the occurrence of electric field threshold-type effects 
and drastic redistribution of proton charge. Electric field applied to correlated chains induces first 
a formation of ordered dipole structures for lower field strength, and than, with a further field 
strength increase, a stabilization of states with Bjerrum D-defects. We analyze the main factors 
responsible for the formation/annihilation of Bjerrum defects showing the strong influence of the 
complex interplay between reorientation energy, electric field and temperature in the dynamics of 
proton wire. 
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I. INTRODUCTION 



Translocation of protons over long distances has a key importance for biological and 
chemical systems. It is believed that proton migration along the chains of water molecules 
formed between the interior of proteins and the solvent, establishes electrochemical potential 
gradients playing an important functional role^^. Experimental evidence indicates that the 
dominant mechanism responsible for proton transport in transmembrane proteins (for in- 
stance, in bacteriorhodopsin of Halobacterium halohiuw^ and in gramicidin A channels^"^) 
is the diffusion of H"^ ions which is faster than the hydrodynamic flow of hydronium species 
(HaO)^. Especially at low hydrogen concentrations in channels, proton conduction is deter- 
mined by a two-stage Grotthuss-type mechanisnt^i^ shown schematically in Fig. 1(a). The 
first stage involves the intrabond proton tunnelling along the hydrogen bridge which is con- 
nected with the formation and transfer of ionic positive (HsO'^) and negative (OH^) charged 
defects. To sustain a fiux of H"*" in such proton wire, the inter-molecular proton transfer due 
to the reorientations of molecular group with proton is assumed. Reorientation motion leads 
to the breaking of the hydrogen bond (so-called orientational Bjerrum L-defect) and location 
of proton between another pair of molecular groups^. Consequently, the reorientation step 
in the presence of the second proton may induce high-energy configuration with both of the 
protons shared by two adjacent oxygen ions (Bjerrum D-defect). 

Unlike the translocation of monovalent ions like Cs^, Na"^ or K"*" via gramicidin requiring 
the net diffusion of the whole water column in the channel, the existence of the Grotthuss- 
type selective migration of through the H-bonded chain is supported by the absence of 
streaming potentials during H"^ permeation^'^. In contrast to the bulk water, the reorienta- 
tion motion in one- dimensional water wire involving a migration of Bjerrum defects with the 
period of reorientations about 10^^° s, is much slower than the proton intra-bond hopping 
(~ 10~^^ s)i2*ii. This is also closely related to the fact that the mobility of Bjerrum defects 
(~ 10"'^ cm^V~^s~^) is much lower than that of the ionic defects (~ 10~^ cm^V~^s~^)^. 
Moreover, as appears from the results of molecular dynamics simulationai^, the translocation 
of the ionic defects in proton wires is almost activationless process, whereas the orientation 
defects involving an activation energy of about 5 kcal/mol in the chain containing eight water 
molecules, constitute a limiting step for the proton migration. Besides the orientation de- 
fects, the recent experiments indicate that another possible rate-limiting step for the proton 
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migration in gramicidin channels can be at the membrane- channel /solution interfaceii*^^. 

As was pointed in^, the experimental analysis of the proton flow in bioenergetic proteins 
and the mechanism of proton translocation is very difficult because of its intrinsically tran- 
sient nature. To shed more light on the microscopic nature of the proton transport and 
to analyze the influence of quantum effects and interaction with proton surrounding, the- 
oretical modeling remains to be essentially important. Recently, much attention has been 
focused on the theoretical studies of the dynamics of ionic defects using soliton model a^^i^^ 
and molecular dynamics simulation3iiiiS*i2iS£. Proton transfer in water was shown to be 
strongly coupled with the dynamics of local environment, and the density of ionic defects 
was found to increase exponentially with the increasing temperature^. However, since the 
concentration of slow Bjerrum defects in water solutions is much higher [cb = 2 ■ 10~^ at 
— 10° C) than that of fast ionic (c/ = 3 ■ 10~^^)^, the investigations of the reorientation 
step of proton migration are necessary for the better understanding of the proton transport 
process. 

The goal of the present work is to study proton wire containing a flnite number of wa- 
ter molecules by the use of quantum statistical mechanics methods which are extremely 
effective in the description of the collective nature of the proton transfer and in the quan- 
tum treatment of the light H nuclei^^. To describe correctly the proton transport process, 
we employ here two-stage proton transport model^^ incorporating quantum effects such as 
proton tunneling and zero-point vibration energy. In earlier paper a^^i^^i^^ we applied the 
two-stage model to analyze the effect of coupling between protons and molecular group vi- 
brations on proton conductivity in infinite H-bonded chains and proton-conducting planes. 
In particular, it was shown that the proton-lattice vibration interactions induce structural 
instabilities and charge ordering in system^^, whereas the Grotthuss-type transport mech- 
anism manifests itself in nontrivial temperature- and frequency dependences of the proton 
conductivity2ii2Si2&. 

In this work we analyze the influence of proton-proton correlations, comparing two dif- 
ferent protonated chains containing one and two excess protons respectively. We flnd that 
the reaction of protonation of water chain is extremely sensitive to the reorientation energy 
barrier of proton motion and the barrier for the chain protonation. We show that the in- 
crease of the reorientation energy results in the drastic decrease of the proton charge density 
at the boundary between the chain and surrounding with consequent localization of protons 
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near the inner water molecules. As appears from our modeling, the application of external 
electric field induces the step-like threshold-type effects with the ordering of proton charge 
and stabilization of Bjerrum D-defects in the wire. We analyze the temperature dependence 
of proton polarization and D-defect concentration, and examine the role of the interplay 
between different factors (such as orientation energy, external field and temperature) in the 
dynamics of Bjerrum defects. 

II. MODEL SPECIFICATION AND METHOD OF CALCULATION 

To model a proton wire, we consider a linear chain containing hydrogen bonds and 
/ = 1, 2, . . . , A^, A^+1 molecular groups. The outer left (/ = 1) and right (/ = A^+1) molecular 
complexes mimic the surrounding of the proton wire and differ from the inner {1 = 2,..., N) 
water molecules. The transport of an excess proton through the wire proceeds via the 
following two steps: 

(i) proton can be transferred within a H-bond (process shown by short arrows in Fig. 1(a)) 
which is modelled by a simple double-well potential, with the corresponding energy barrier 
Qt for the proton transfer between the two minima: 

N 

Ht = ^{clcib + 4^Cia), (1) 
1=1 

where c^(Qiy) are the operators of the proton creation(annihilation) in the position (/, u) 
(the index u = {a, b} denotes the left/right position for the proton within the H-bond); 

(ii) a water molecule together with covalently bonded hydrogen ion can be rotated, and this 
process causes the breaking of the H-bond and location of H"*" between two another nearest 
water molecules of the wire (process depicted by long arrows in Fig. 1(a)): 

N 

Hr = nR J2i4+l,aClb + 4CZ+I,a), (2) 

1=1 

where Qr is the effective energy barrier for the proton hopping between the states \l + l,a) 
and \l,b) (reorientation of the Z-th molecular group together with proton). As is shown in^I by 
the computation of the proton mean-force potential, this transition between donor-acceptor 
and acceptor-donor states reverses the chain dipole moment and requires a substantial energy 
barrier about 5.5 kcal/mol for the whole chain containing up to eight water molecules. 



4 



Besides the transport process, we incorporate the following two types of interactions 
between protons in the chain: 

(iii) different short-range configurations of the protons near an inner water molecule as well 
as an outer molecular group can appear due to the different nature of bonding (shorter 
covalent or longer H-bond). The energies of possible configurations (shown in Fig. 1(b)) are 
described by the following terms: 

Hi = ei{l - nia) + wiTiia, Hn+i = eN+i{l - nm) + WN+iUm, (3) 
Hi = w'{l - n;+i,„)(l - nib) + wni+i^aUib + e{l - nib)ni+i^a + enib{l - rii+i^a)- 

The parts Hi and H^+i describe the energies of the boundary proton configurations near 
the left and the right surrounding molecular groups (we assume for simplicity ii = ejy+i = e 
and Wi = wn+i = w in the boundary configurations shown in the upper scheme of Fig. 1(b)). 
The terms Hi {1 = 2,..., N) contain the configuration energies for the water molecules in 
the interior of the wire (the lower part in Fig. 1(b)). Here the proton occupancy operators 
nil, = c^ciu = {0, 1}; 

(iv) a strong repulsion between two nearest protons shared by two neighboring oxygens (so 
called Bjerrum D-defect) with a repulsion energy U is represented by the term: 

N 

Hc = U^nianib. (4) 

1=1 

In our following analysis we use the value of U ^10 kcal/mol corresponding to the energy 
of relaxed D-defect estimated in2& on the basis of quantum chemical calculations. 

To model a field exerted by the surrounding, we apply an external electric field of a 
strength E to the proton wire, which is described by the following term 

N 

-CpE ^ Ri^ni^, (5) 

1=1 

f — {a.b} 

where Rii, is the coordinate of the proton position (/, u) with respect to the center of the 
chain, and Cp denotes the proton charge. 

In order to analyze the dynamics of the proton wire embedded in the surrounding under 
the influence of the field, as well as the effect of rotational motions of covalent groups with 
proton, we will focus our attention on the polarization of the proton wire defined here as 

N 

P = ep ^ Riu{niu), (6) 

t' — {a,6} 
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where (. . .) denotes the statistical average with respect to the system energy ((TJIS)). The 
average probabihties {riiy) of occupation of the position (Z, v) by proton describe the distri- 
bution of the proton charge in the wire, and thus is another very important characteristics 
to track the proton migration. 

To calculate exactly the above-mentioned statistical averages, we need to know the quan- 
tum energy levels determined by the energy (PEI). This can be done by a mapping of the 
proton states {l,v) on the multi-site basis \i) = \nia,nib, . . . ,nNb)- Then, using the projec- 
tion operators X**' = acting on the new basis \i) we rewrite the system energy 
in a convenient form (see Appendix): 

H = H'^®H^®---® H^^. (7) 

Each term H'^^ in ((7j) corresponds exactly to Hp protons in the chain [rip = "^^i ui^ii^) ~ 
0, 1, . . . , 2N). This means in fact that the mapping on the states \i) allows to decompose 
the terms (PIH)) and analyze the cases of different number of protons in the wire separately. 

The energy barrier for the protonation of the chain is described by the parameter A 
which appears in (|7j) after the decomposition procedure (see Appendix). As follows from 
the definition ()A8|) . A is the difference between the energies of the proton attraction to the 
boundary and to the inner water molecules. As the PMF-studies of the H-bonded chain 
dynamica^l show that the inner H-bonds are stronger (shorter 0-0 separation distances) 
than the outer H-bonds, it is reasonable to consider below the case A < (we take |A| = 
0.85 kcal/mol in our numerical calculations), when the proton is attracted to the surrounding 
and needs to overcome the boundary energy barrier |A| to protonate the water chain. 

The parameter J = w + w' — 2e (see Fig. 1(b)) is related to the effective short-range 
interactions between the protons near the water molecule. It describes, in fact, the energy 
of the formation of the pair of ionic defects: 1+=H30^ {w — e) and I_=OH~ {e — w') from 
two water molecules at the dissociation reaction (2H2O H^O^+OB.^). Since the value of 
J is about 22 kcal/moii^i^ and is more than twice as much as U, we exclude in our following 
analysis an appearance of the pair of ionic defects in the system. 

Due to the two types of motions we have two different contributions to the proton 
dipole moment: the orientational part fir = CpRr and the transfer part fiab = GpRab where 
Rah = Rib — Ria denotes the distance H-H between the two nearest proton positions of the 
double-well H-bond. In our calculations, we use the values fiab = 3.5 D and fiab = 4.5 D 
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corresponding to the moderately strong H-bond with an 0-0 distance Rab + 2/?^ = 2.6A 
and the covalent 0-H bond of a length = 0.94 A. 



III. CASE STUDY: ONE EXCESS PROTON IN CHAIN 

As the starting point, in this section we mimic the situation when one proton is moved 
towards the water chain embedded into the solvent. To examine the behavior of the proto- 
nated chain with N H-bonds and n = 1 excess proton, we consider part of the energy 
given by ()A4j) . Since the zero-point vibration energy for protonated chains is larger than the 
potential energy barrier for the proton transfer between two shared oxygens^, the quantum 
tunneling is not required for the intra-bond transfer. Thus, in our modelling we set 
Qj. ^ 0. With this assumption, the energy levels of can be found exactly: 

Al,2 = A±(^{N- l)flr + y/iafc^ E, 

. ±i{^ + ^ir)E±p (i = 1,3,..., AT -2), forodd 
K-,2N = { (8) 
±i{tY + l^r)E ±p (i = 0,2,...,A^-2), for even A^ 



where p = ■>/ /i^-E^ + 

To analyze the role of Qpt we consider first the case without external field {E = 0). 
Depending on the value of Qr, two different regimes may be stabilized in the system. In 
the first small-Qji regime, the two lowest energy levels Ai_2 correspond to the superposition 
of the two boundary states 

I10...00) and |00...01) (9) 

with the proton located in the surrounding near the left or the right outer molecular groups. 
In the second large-Qn regime the proton is shared between the inner water molecules of the 
chain and the ground state of the system corresponds to the superposition of the states 

|010...00),|0010...00) I00...010) (10) 

with the energies A4, Ag, . . . , A2Ar. The "critical" value = —A separating these two 
regimes, reflects the transition of the proton from the surrounding to the states where the 
proton is shared by the chain water molecules, which corresponds to the protonation chemical 
reaction. Fig. 2 shows the variation of the average occupancies of proton sites with Qji for the 



chains containing N = 2 and = 3 H-bonds. For low temperatures (see the case T = 30K) 
the boundary proton occupancies (riia) = {riNb) drop drastically to zero at Qr = whereas 
the occupation numbers of the central positions (nib) = ■ ■ ■ = {riNa) increase up to the value 
2(Af"-i) ? reflecting the redistribution of collectivized proton between the inner sites in the 
wire. It should be noted here that Q*^ reflects the change of the ground state of the system 
and is determined as the solution of the equation Ai^2 = -^4^6, ■■ ,2N at = which does not 
depend on temperature. However, as all statistical averages, the average proton occupancies 
(for example, of the states = :^(|l5&) + |2,a)) and |2,a) = -^(^11)^) + |2,a)) where 
1 1,6) and |2, a) are the diagonalized states corresponding to A4 and A3 respectively), are 
temperature dependent. Thus the value of ^^^(T) ^ —A + A;Tln2 where (riia) = {n2b), 
for T 7^ is not equal to (see Fig. 2, case T = 300K). This difference shows that 
the inner proton states (fTUI) are stabilized already at lower Qji = < Q*j^{T), although 
the occupancies of the inner positions a.t = Q*^ are still slightly lower then of the outer 
due to the temperature- induced fluctuations. As T 0, the fluctuations decrease and 

The effect of the proton localization inside the chain is supported by the results reported 
iii2L2i showing that in H-bonded finite chains, the excess charge is best solvated by the 
central H-bonds. However, as results from the presented above analysis, the effect of proton 
localization is drastically influenced by the competition between two different tendencies: (i) 
for small Qr, the proton is located near the surrounding/ wire interface due to the nonzero 
protonation barrier A; (ii) to overcome the barrier between the surrounding and the wire, 
the reorientation energy should be sufficiently large {Qr > —A) in order to stabilize the 
inner proton configurations. These conclusions show that in general, these two different 
factors (interface barrier and orientations) can be rate-limiting for the charge translocation 
and proton conductivity of the wire. As was shown in-, the effective reorientation barrier 
can be influenced by temperature factor or applied voltage (for example, the reorientation 
rate of the wire decreases exponentially with T decrease). Thus, one can also expect that 
the increase of the temperature can result in the lower orientation barrier Qr, delocalization 
of proton and consequently in higher values for the proton conductivity through the channel. 
However, more detailed theoretical analysis is needed to understand better the role of the 
interface in the behavior of the conductivity. 

We turn now to an analysis of the proton translocation directed by the external field (the 
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case 7^ 0). Fig. 3 shows the field-dependences of P and (rii^) for the chain containing 
N = 2 H-bonds. We note that the behavior in the first small- and in the second large-fi^ 
regime differs drastically. In the small-fi/j regime the polarization increases smoothly with 
E approaching finally its maximal saturation value Pmax = fJ-r + A^afe (Fig- 3(a), inset). In 
contrast to this, in the large-Q^ regime the field dependence is rather nontrivial: first, for 
E < Ethresh, P changes very slowly, and than, at ~ Ethresh a strong increase of P to 
P = Pmax is observed in Fig. 3(a). This rapid step-like change of the proton polarization 
refiects the threshold-type effect where the threshold electric field value at low T is given by 

P _ A ■ /in + VA^/i^ + ^Kf^H - f^r) ( - , \ 

^thresh — o 5 5 [f^H — l^r + l^ab) 

and does not depend on the chain size {Ethresh ~ 0.15-10^'' V/cm for the chains with = 
2 as can be observed in Fig. 3(a)). As we see in Fig. 3(b), the proton charge translocation 
under the infiuence of the field in this case proceeds not smoothly, but has a step-like 
character. As results from (fTT|) . the threshold value Ethresh (which is needed to overcome a 
barrier for pumping between the inner localized states ()10|) and the boundary state |00 . . . 01) 
in the direction of field) increases for larger Qji (Fig. 3(a), the cases with = 1.5 kcal/mol 
and Qr = 2.15 kcal/mol). This implies that the conductivity of protonated chains can 
drop with an increasing which can occur in system for example due to the temperature- 
induced fluctuations of Qr. However, as was shown ini2^, the reorientation rate increases 
at increased voltage, which corresponds in our case to the smaller values of Qr for E ^ 0. 
Thus, we can expect that the external field-induced lowering of the orientation barrier for 
the proton translocation results in the increase of the proton conductivity in the wire. 

The drastic change of P at E = Ethresh leads to a strong qualitative difference in the 
temperature shapes of the polarization profiles shown in Fig. 4(a) for the large-Qn regime. 
For E < Ethresh and Qr > Q*^, the excess proton is located in the central sites, and the 
increase of T results in a disorder-induced transfer of the proton from the inner positions 
to the chain boundary giving the increase of P at 300 K as compared to 50 K (Fig. 4(a), 
cases E = 0.02 ■ 10"^ V/cm and E = 0.1 ■ 10"^ V/cm). As the proton is located in the outer 
state |00. . .01) for E > Ethresh (corresponding to P = 1 for T — ^ 0), the dominant effect 
of T in this case is the the disorder-induced proton redistribution between all sites in the 
chain leading to the lowering of total polarization (Fig. 4(a), cases E = 0.2 ■ 10^ V/cm and 
E = 0.3 ■ 10^ V/cm). The profiles of the polarization for Qr below fi^j shown in Fig. 4(b) 
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appear to be very similar to the high-field profiles in Fig. 4(a). Since for < Q*^ the proton 
is located in the outer states near the chain boundary already at low E, the increase of 
T suppresses the polarization due to the increasing proton disorder. 



IV. ROLE OF PROTON-PROTON CORRELATIONS 



In order to examine the influence of proton correlations, we consider next the translocation 
of two excess protons in the wire which is described by the part of the total energy ()A4j) . 
Since the presence of two protons in wire may lead to the formation of Bjerrum D-defect, 
the energy for the chain with N = 2 H-bonds given by the expression ()A6|) . contains the 
terms with the energy U of the repulsion between two nearest-neighboring protons. 

Analogously to the 1-proton wire, we analyze first the behavior of the chain without the 
electric field. For N = 2 and Qt ~ the energy levels found from ()A6jl have the following 
form: 



Ai,2 = -{U ±q-)+ hhE, 



A3 = 
^5,6 



A, A. 



-A + J, 



(12) 



1 



{U±q+)-fiHE, 



(9. 



(^7 T 'i^rEf + An\) 



and correspond to the following states of the wire: 

P- 



1) 

2) 

3) 
5) 

6) 



:|6) + 
:|6)- 



p- 



|8), |4) = |9), 

HM\ , P+ 



|7), 
|7), 



(13) 



p+ 



m + 



1 10) 



p^ 



111), 



n 



R 



111) 



with p± = ^{U ± q±) =F fJ'rE. 

To study the influence of Qr, we analyze (|T^ and (fT!?jl for = assuming J U 
and neglecting in this way by the formation of ionic defects described by the configuration 
|9) = |0110). Similarly to the 1-proton wire, the two different regimes can exist in the 
system depending on the value of the reorientation energy. In the first small-VL^ regime 
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(for ^Iji < ^*^^^ = VA^ — UA), each proton is located near the outer molecular group and 
the state |8) = |1001) has the lowest energy A3 = A. As increases and approaches 
the "critical" value ^*^^\ the transition to the large-Qji regime occurs. In this regime (for 
> n]^^^) the lowest energy levels A2 = Ae = |(?7 - q) (q = q+{E = 0) = q_{E = 0)) 
correspond to the states |2) and |6) in (fT!^ . with one proton located in the interior of the 
wire. However, in contrast to section 11111 the transition between these two regimes is U- 
dependent, because fi^^'^ contains the energy of the D-defect U. Fig. 5(a) shows the variation 
of the proton occupation numbers {rii^) with Qf> for U = 9.4 kcal/mol and U = 1.9 kcal/mol 
(plotted in the inset). The "critical" value ^*r^^ ~ 3 kcal/mol for the repulsion energy 
U = 9.4 kcal/mol is larger as compared with fi^^'"* ~ 1.5 kcal/mol for U = 1.9 kcal/mol. So 
far as the repulsion energy U is significant, ^l*^^^ > fijj. However, as f/ — > 0, ^*r^^ approaches 
the "critical" value Q*^ for the one-proton case . As we can see from ()13|) and ()A1|) . the 
ground states |2) and |6) in the large-fi/j regime are represented by the superpositions of the 
normal configurations |7) and |10) and the states |6) and |11) containing the D-defect. Thus 
the transition to the large-fi/j regime stabilize D-defects inside the chain. The formation 
of the D-defect states is clearly observed in Fig. 5(a) where the occupation numbers of the 
D-defect states |6) and |11) significantly increase for > nf\ 

Analogously to the one-proton case, the temperature fluctuations lead to the slight 
temperature-induced increase of the value ^1*^^\t) (corresponding to = = ng), as com- 
pared to fl*^"^^ where the states |2) and |6) are already stabilized. In fact, for Qr > Q*^^\t) 
the states corresponding to protonated chains with significant concentration of D-defects 
prevail {112 = tiq > n^), whereas the states with the protons located at the boundaries in the 
surrounding dominate for Qr < Q*^^\t) {n^ > = rig). For small VLr <^ f/, VL*^^ can be 
given by 

0*J2) r I 8A;Tf/2([//2- A)] ^^^^ 



0*(2) 



A2 - f/A 

Since the line Vt^ (T) found from (fT^ . is tilted with respect to T in the state diagrams 
(T, Vt^\T)) (Fig. 5(b)), the effect of temperature for the chain in the large-fi/j regime 
is crucial: with the increasing T, the temperature fluctuations can destroy the D-defects 
and redistribute the proton charge between the other chain sites. See for example the case 
of Vt*^\T) = 3 kcal/mol and f/ = 9.4 kcal/mol plotted in Fig. 5(b) where the D-defects 
annihilate at T ^ 100 K. 
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As U oo, the weight constants for the D-defect states in (fT^ become smaller: 
. ^ = . , = 0. Thus, the contribution of the D-defect-states to the sta- 
ble wire configuration goes down as ^ for the stronger proton repulsion U (see for the 
comparison {uq) = {rin) for different U plotted in Fig. 5(a)). 

The fact that the variation of temperature can lead to formation or annihilation of the D- 
defects is also observed in the T-dependence of the proton polarization. Note that especially 
for weak external field E, the behavior of P{T) in the small-f2;^ (Fig- 6(a), E = OA-lO'^V/cm) 
and in the large-f2/j regime (Fig. 6(b), E = 0.1 ■ lO^V/cm) is drastically different. In the 
first case, at low T, the predominantly occupied symmetric ground state |10 . . . 01) has the 
total polarization P = 0. However, with T increasing, protons tend to occupy the excited 
states with non-symmetric charge distribution that results in an increase of P. Fig. 7(a) 
demonstrates that the population of all excited states, in particular those containing D- 
defects (Fig. 7(a), inset), grows with T. Although the concentration of the D-defect states 
|0011) and |1100) is of 3-4 orders lower than that of the normal states (see Fig. 7(a), inset), 
it significantly increases up to 1-2 orders with the temperature increase from 50 to 300 K. 

In contrast to this, in the large- 17 /j regime the protons and stable D-defects migrate in the 
direction of applied field ii^ for T ^ giving a non-zero P (Fig. 6(b)). As T increases, the 
population of the excited non-defect state with the protons redistributed at the boundaries 
grows (see Fig. 7(b)) which gives the lower chain polarization. The D-defects, located near 
the end of the chain for finite E (|0011)), are redistributed between another chain positions 
with T, which is observed in Fig. 7(b) showing a decrease of |0011)- together with a slight 
increase of the |1100)-state population at T = 300 K as compared to lower temperatures. 

We study now the electric field effect in correlated chains. Fig. 8 shows the variation of 
polarization and redistribution of protons with increasing E. Consider first the small-fii^ 
regime. In distinct to the 1-proton wire, where the polarization increases smoothly to its 
maximal value Pmax (Fig. 3(a), inset), we observe here two different threshold effects. The 
first transition from the state |8) of ()AH) (the ground state of the wire in the small- regime 
aX E = Q) to the state |10) (where both of the protons are ordered in the right position of 
each H-bond in the direction of the field) occurs at the threshold field value 

" (A^-l)(2/i, + /i,fe)-2/i/ ^^^^ 
The distribution of the occupation probabilities (rig) and (riio) for the states |8) and |10) 
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is plotted in Fig. 8(b). We observe ai E = Ei the abrupt increase of {uiq), while at the 
same field value (ng) drops to zero. Furthermore, we conclude from (fT3j) that the value Ei 
lowers with the number of the water molecules in the chain. This effect can be observed 
in Fig. 9 where the jumps of the polarization are plotted for different A^. Finally, for very 
long water chain (A^ — oo) Ei 0. In contrast to the strong A^-dependence of Ei, the 
second threshold effect appears at 

^= = I; (I'') 

essentially due to the proton correlations and does not depend on the chain length. The 
strong increase of P at E = E2 shown in Fig. 8(a) and Fig. 9 is related to the second 
drastic redistribution of the proton charge in the wire. As can be observed in Fig. 8(b), at 
E = E2 the occupation probability (^11) of the D-defect-state |11) drastically increases to 1, 
whereas (riio) drops to zero. Thus, as resulted from our model, the formation of D-defect in 
external electric field has a step-like character proceeding via the threshold mechanism. In 
the laige-Qpt regime, where the protons are stabilized at the inner water molecules already at 
E = 0, the first threshold phenomenon aX E = Ei, observed for the small-fi/j case, does not 
occur. However, the transition at E = E2 with the increase of the D-defect concentration 
appears in this regime similarly to the regime of small Qr, that can be observed in the P- 
profile for = 3 kcal/mol shown in Fig. 8(a). Note that the effect of the increasing double 
occupancy due to membrane potentials has been observed in the current / concentration plots 
in gramicidin channels^, thus supporting our main conclusions about the role of the external 
electric field. 

The discussed above formation of the D-defects for > Q*^^ in the high electric field 
results in the increase of P for lower temperatures as shown in Fig. 6(b). Basically, the 
essential effect of T observed in the P(T)-profiles in Fig. 6, is the suppression of the total 
polarization due to proton disorder. However, the shapes of the polarization in Fig. 6(a) 
are drastically different for E < Ei and E > E2. For low fields {E < Ei) the polarization 
first increases (refiecting the fiuctuation-induced expansion of proton charge from the outer 
symmetric positions |10. . .01) with P = to the inner positions of the chain accompanied 
by the formation/annihilation of D-defects), and then smoothly decreases due to the disorder 
effect. In contrast to this, as the increasing electric field induces the step-like formation of 
D-defects in the small-fij^ regime, the temperature behavior of P in this case is similar to 
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the the large-Qn case (compare Fig. 6(a) and Fig. 6(b) with E = 1 ■ 10^ V/cm) showing the 
smooth disorder-induced decrease of P with T. 

We also note that the stable configurations with double proton occupancy require the 
additional reorientation steps for the proton translocation and can result in the smaller 
values for the proton conductivity. This fact has been observed in the measurements of 
the proton conductance in two different stereoisomers of the gramicidin^, thus supporting a 
possibility of stabilization of the D-defect states in proton wires. 

V. SUMMARY 

In this work we studied the process of proton translocation in ID-chains mimicking pro- 
tonated water channels embedded in surrounding. We have analyzed the role of the reorien- 
tation motion of protons, as well as the effect of electric field and proton correlations on the 
chain dynamics. We have shown that the increase of the reorientation energy results in the 
transition to the large-fi^j regime characterized by the transfer of the proton charge from 
the surrounding to the inner water molecules in the chain. The process of proton migration 
along the chain in the external electric field has the step-like character leading to the ap- 
pearance of the electric field threshold-type phenomena with drastic redistribution of proton 
charge. The correlations between protons in the chain increase the "critical" reorientation 
energy necessary for the transition into the large-fi/j regime, where the protonated chain 
contains a finite concentration of Bjerrum defects. The temperature fiuctuations induce a 
slight increase of fl*j^{T) separating the state with the protons located in surrounding near 
the outer groups, and the protonated state with D-defects. For the correlated chains, this 
temperature dependence of the "critical" reorientation energy can lead to the redistribution 
of proton charge and annihilation of D-defects with increasing T. The electric field applied 
to the correlated chains induces first the formation of ordered dipole structures for the lower 
E values, and than, with the further E increase, the stabilization of the states with the 
Bjerrum D-defects. 

Generally, the increase of temperature suppresses the total polarization in the chain due 
to the increasing disorder. However, especially in the low electric fields, the shapes of the 
temperature profiles of the polarization appear to be drastically different in the small- and 
large-f2ij regimes demonstrating the complex interplay between the reorientation energy and 



14 



temperature. 

Finally, as follows from our analysis, the following factors strongly influence the formation 
of Bjerrum defects: (i) the high electric fields can form the defects and pump them in the 
chain in the direction of field; (ii) the increase of the orientational energy barrier leads to the 
stabilization of D-defects; (iii) the increase of temperature in the large-f2^ regime results in 
the formation/annihilation of D-defects, whereas for small Qr the concentration of D-defects 
significantly increases up to 1-2 orders at the room temperatures as compared to the low 
T 50 K. 
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APPENDIX A: DECOMPOSITION OF THE PROTON STATES IN THE SYS- 
TEM WITH N = 2 H-BONDS 

We demonstrate below the procedure of the mapping in the system with N = 2 H-bonds 
on the multi-site states. For N = 2 the basis \i) includes 2^^ = 16 states \nia,nib,n2a,n2b)'- 

|1) = |0000), |2) = |1000), |3) = |0100), 

|4) = |0010), |5) = |0001), |6) = |1100), (Al) 

|7) = |1010), |8) = |1001), |9) = |0110), 

|10) = |0101), 111) = |0011), . . . |16) = |1111). 

We can derive the relations between ci^^ and X"' = 

ci,u = 5^(^|Q,,|J)X^^ (A2) 

where the expectation numbers {i\ci^u\j) can be found using the usual antisymmetric rules 
for Fermi-operatora^S. Specifically, for the case N = 2 the expressions ()A2|) yield: 



Co,a 


= Xl'2 


+ X3'6 


+ X^'^ + X^'S + X9'12 + Xl0'13 ^ ^11,14 


^^15,16 




Cl,a 


= X''^ 




_ ji^3,9 _|_ Ji^5,ll _|_ Ji^6,12 _ Ji^8,14 _ j|^10,15 


^^13,16 


(A3) 


C0,b 






+ X^'9 + X^'IO - X^'12 _ ^8,13 ^ ^11,15 


_ j|^14,16 




Cl,b 


= Xi'5 


-X2'8 


- X^'io _ x^'ii + X6'i3 + x^'i^ + X^'is 







Using the relations ()A3|) and the fact that X**'X"' = ^iqX^^' (due to the orthogonality of 
the states |z)), we decompose (013) in terms of X"' operators into the following 5 terms: 

H = Hl® Hi ® Hi ® Hi © Hi, (A4) 

where 

— ^0 

-"2 — 0,25 

Hi = {A + {flr + flab)E)X^'^ +firE{X'''^ -X^'^) + 

(A - (/i, + fiab)E)X'>'' + QriX''^ + + (A5) 
fir(X^'' + X^'^) + nii{X^'^ + X^'3) + al, 

16 



nniX^^' + + + + aj, (A6) 

= {U + J-A + {fIr + flab)E)X^^''^ + {U + firE)X^^'^'' + 

{U - firE)X^^''^ + {U+J-A-{fIr + fiab)E)X^''''^ + (A7) 
fiT(Xl2,13 ^ ^13,12) ^ ^^(X^^'l^ + X^S'l^) + 

Since the parameter 

A = {w-e)-{e- w') (A8) 

in ()A5|) - ()A7|1 is the difference between proton configuration energies at the boundary = 1 
or Z = 1 surrounding molecular groups), and at the inner (2 < / < A^) water molecule, it 
describes, in fact, the energy barrier for the protonation of the water chain. For our analysis 
A has the key importance, because the other energy constants in ()A5|) - ()A7p 

a2 = al — {w — e) — (e — w'), al = al — (w — e), 

al = al + {w — e), a2 = al + {w — e) + {e — w'), 
al = w + e + e 

which appear due to the boundary effects, are independent of the proton location in the 
wire and thus do not influence the statistical characteristics like ©• 

By the similar way the energy can be rewritten for the systems with any finite value 
of N. 



^ Gennis R B 1989 Biomembranes: Molecular Structure and Functions (New York: Springer- 
Verlag) p 235 

^ Lanyi J K and Pohorille A 2001 Trends in Biotechnology 19 140 

^ Stoeckenius W, Lozier R H, and Bogomolni R A 1979 Biochim. Biophys. Acta 505 215 

17 



^ Nagle J F and Mille M 1981 J. Chem. Phys. 74 1367 

^ Levitt D G, Elias S R, and Hautman J M 1978 Biochim. Biophys. Acta 512 436 

6 Cukierman S 2000 Biophys. J. 78 1825 

^ Chen M S, Onsagcr L, Bonner J and Nagle J 1974 J. Chem. Phys. 60 405 

^ Bjerrum N 1951 K. danske Vidensk. Selsk. Skr. 27 1 

^ Akeson M and Deamer D W 1991 Biophys. J. 60 101 

1° Nagle J F, Mille M, and Morowitz H J 1980 J. Chem. Phys. 72 3959 

1^ Chernyshev A and Cukierman S 2002 Biophys. J. 82 182 

Eisenberg D and Kauzmann W 1969 The Structure and Properties of Water (Oxford: Oxford 
University Press) 

13 de Godoy C M G and Cukierman S 2001 Biophys. J. 81 1430 

14 Pomes R and Roux B 1998 Biophys. J. 75 33 

15 Savin A V and Zolotaryuk A V 1991 Phys. Rev. A 44 8167 

1^ Zolotaryuk A V, Savin A V, and Economou E N 1994 Phys. Rev. Lett. 73 2871 

1''' Tuckerman M, Laasonen K, Sprik M, and Parrinello M 1995 J. Chem. Phys. 103 150 

1^ Marx D, Tuckerman M, Hutter J, and Parrinello M, 1999 Nature 397 601 

1^ Schumaker M F, Pomes R and Roux B 2000 Biophys. J. 79 2840 

20 Schumaker M F, Pomes R and Roux B 2001 Biophys. J. 80 12 

21 Davydov A S 1982 Biology and Quantum Mechanics (Oxford: Pergamon Press) p 25 

22 Stasyuk I V, Ivankiv O L and Pavlenko N 1997 J. Phys. Studies 1 418 

23 Pavlenko N 2000 Phys. Rev. B 61 4988 

24 Pavlenko N 2000 J. Chem. Phys. 112 8637 

25 Pavlenko N 2000 Phys. Stat. Sol. (b) 218 295 

26 Pavlenko N I and Stasyuk I 2001 J. Phys.: Condensed Matter 13 4081 

27 Pomes R 1996 J. Phys. Chem. 1996 100 2519 

28 Hassan R, Campbell E 1992 J. Chem. Phys. 97 4326 

2^ Davydov A S 1965 Quantum Mechanics (Oxford: Pergamon Press) 

30 Rao CNR 1972 Theory of Hydrogen Bonding in Water, in: Water: a Comprehensive Treatise 
vol 1 (New York: Plenum Press) 

31 Pomes R 1999 Isr. J. Chem. 39 387 

32 Phillips L R, Cole C D, Hendershot R J, Gotten M, Cross T A, and Busath D D 1999 Biophys. J. 

18 



77 2492 



19 






/ 

• 


\ 


\ 




w 




w 


A. 




A, 


A. 


w' 


w 


s 


s 



(b) 



FIG. 1: (a) Schematic presentation of proton wire, arrows indicate a possible path of proton 

migration along the chain. Full circles denote water molecules and open circles are the possible 

positions for excess proton, (b) Scheme of possible proton configurations near the outer surrounding 

groups (the upper part) and the inner water molecules of the wire (the lower part). 
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FIG. 2: Proton site occupancies for different reorientation energies in the chain with two 
hydrogen bonds containing one excess proton. The inset shows the redistribution of proton charge 
in the chain with iV = 3 H-bonds. 
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FIG. 3: (a) Proton polaxization vs electric field in the small-Ofl (shown in the inset)- and large- 
regimes for T = 30K, and (b) average site occupancies vs E in the chain containing N = 2 
hydrogen bonds and n = 1 proton for CIr = 1.5 kcal/mol and for different temperatures. 
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FIG. 4: Proton polarization vs temperature (a) in the laige-^u regime and (b) in the smah-OR 
regimes for different values of applied electric field E. 



23 




3.5n 

3.4- 

_ 3.3- 

o 

E 

^ 3.2- 
o 



a 



3.1 - 
3.0- 
2.9- 



1.80 
-1.75 
E 1.70 
8 1-65 

J. 60 



1.55 
1.50 



U=1.9kcal/mol 

protonated chain 
with D-defects^-^ 



protonated chain 
with D-defects 



protons at the boundaries 

50 100 150 200 250 300 



(b) 



protons at the boundaries 

50 100 150 200 250 300 

T, K u=9.4 kcal/mol 



FIG. 5: (a) Proton charge distribution vs Qr in the H-bonded chain {N = 2) containing two 
protons (n = 2) for T = WOK. The inset shows the variation of proton charge with J7r for 
U = 1.9 kcal/mol. (b) State diagrams (T, fi^) for N = 2, n = 2 indicating the regions of stability 
of the protonated chain with D-defect and the states with the protons localized at the boundaries. 
The inset shows the state diagram for lower U = 1.9 kcal/mol. 



24 



Q. 

5: 

c 
g 

N 



E=0.4 10' V/cm 
E=0.5 10'V/cm 
E=1.0 10'V/cm 




50 100 150 200 250 300 350 
(a) T, K n =0.43 kcal/mol 



Q. 

qI 

c 
o 

' 

CO 
N 
'l_ 

o 

Q. 

c 
'co 



0.40 
0.35 
0.30 
0.25 
0.20 
0.15 
0.10 
0.05 
0.00 



E=o.i lo'v/cffr 

E=0.4 10'V/cm 
E=1.0- 10W/cm 



50 100 150 200 250 300 350 



(b) 



T, K 



n„=3.0 kcal/mol 



FIG. 6: Proton polarization in the chain with N = 2 H-bonds and n = 2 protons vs temperature 
(a) in the smaU-riij regime and (b) in the large-fi^j regimes for different values of applied electric 
field E. 
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FIG. 7: Proton site occupancies in the chain containing N = 2 H-bonds and n = 2 protons vs 

temperature (a) in the small-r^i? regime and (b) large-f^i? regimes for low electric field E. The 
inset shows the concentration of D-defects vs T in the small-fiij regime. 
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FIG. 8: (a) Proton polarization for different and (b) average site occupancies for ^Ir = 
0.43 kcal/mol vs electric field E in the chain containing N = 2 hydrogen bonds and n = 2 protons 
at T = 30K. 
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FIG. 9: Proton polarization vs E in the chains of different length containing two protons for 
nn = 0.43 kcal/mol and T = 300K. 
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